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Abstract 

We present a semiclassical study of a transport process, the tunneling, in the pres- 
ence of a magnetic field and a dissipative environment. We have found that the problem 
can be mapped onto an effective one-dimensional one, and the tunneling rate is strongly 
affected by the magnetic field, such as a complete suppression by a large parallel mag- 
netic field, an example of the dynamical localization. In such case a small perpendicular 
component of the field, or the dissipation, can enhance the tunneling rate. In the small 
parallel field and finite temperatures the tunneling rate is finite. Explicit expressions 
will be presented in those cases. If viewing the tunneling in the presence of a magnetic 
field as a dissipative tunneling process, by varying the magnetic field and the potential 
one can obtain the dissipative spectral function between the subohmic s = and the 
superohmic s = oo. In combination with a real dissipative spectral function, the effect 
of the magnetic field can map the spectral function from s to 2 — s, with s > 2 mapping 
to s = 0, revealing a dual symmetry between the friction and the Lorentz force. Two 
cases relevant to experiments, the edge state tunneling in a Hall bar and the tunneling 
near the dynamical localization will be discussed in detail. 
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1 Introduction 



Quantum dissipative tunneling has received an extensive study recently]!], 0, because of 
its important role played in physics, chemistry, as well as in biology. In the study of the 
quantum transport in semiconductors, such as in heterojunctions, superlattices, and the 
quantum Hall systems, this problem surfaces again, along with a new feature in addition to 
the usual dissipative effect: the effect of the magnetic field. [0, U This new feature makes 
the tunneling behavior much more rich. The purpose of the present paper is to explore its 
various consequences.^ 

Classically, the equation of motion for a charge particle with mass m and charge q is 



with the potential U, the friction coefficient 77, the magnetic field B, and the speed of light 
c. The fluctuating force is connected to the friction by the fluctuation-dissipation theorem. 
Eq.(l) is the formiliar Langevin equation, and the friction here corresponds to the ohmic 
damping case. We notice here that both the damping and the magnetic field give velocity 
dependent forces which are perpendicular to each other. This suggests a certain symmetry 
between the the friction and the Lorentz force, which can indeed be borne out as will be 
shown below. 

To study the quantum tunneling, the Hamiltonian of the whole system is needed. For 
the Lorentz force, it is straightforward. The corresponding Hamiltonian is 



with the vector potential A determined by V x A = B. The formulation of the dissipation 
within quantum mechanics is not so straightforward. In the present paper the approach and 
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notation of Ref.[l|] to this question will be followed. In that formulation, the dissipation 
is produced by a heat bath consisting of many harmonic oscillators, and the corresponding 
Hamiltonian is 



H 



D 



E 

j;l=x,y,z 



1 2 1 2 



Ri 



Here (R x , R y , R z ) = R 
the spectral function 



(x,y,z). The effect of the dissipative environment is specified by 



(2) 



with I = x,y, z. The introduction of the notation I stands for the case of a possible anisotropic 
dissipative environment. In the present paper, we shall assume the spectral function to have 
the following form 

J l (u) = r ]l u Sl exp{- — \ , (3) 

with uj c the cutoff frequency whenever needed. In accordance with Ref.|IJ, si > 1 is the 
superohmic case, sj = 1 the ohmic case, and < si < 1 the subohmic case. In the ohmic 
damping case, 7]i is the friction coefficient in Eq.(l). 

Now, the total Hamiltonian of our system has the following form: 

\ 2" 
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Pl,j + o m j^j %i 
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(4) 



We further assume the scalar potential U for the charged particle to have the form: 

u = \ m u 2 x x 2 + ]^nuj 2 y y 2 + V(z) , 



(5) 



with V has a metastable point at z — 0. For large enough z, V(z) = —Fz, and there is a 
barrier separating this region and z = 0. The choice of the potential in the x — y plane is to 
mimic the impurity potential. The external magnetic field is tilted: B = (B\\,Q, B±). Here 
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the magnetic field parallel to the x — y plane is taken along the x-direction. In accordance 
with the calculation of tunneling, the vector potential will be taken as 



and it can be shown that results are independent of the choice of gauge, because of the 
periodic boundary condition in the calculation of the tunneling rate. It is clear now we have 
encountered the problem of tunneling in a higher dimensional space: the particle moves in 
3-d, and the effective dimension for the dissipative environment is infinite. The suitable tool 
to deal with such a problem is the path integral method, which will be used below. 

The organization of the rest paper is as follows. In the next section using the path integral 
method we reduce the infinite dimensional problem to an effective one dimensional one, with 
normal and anomalous damping kernels. The situation corresponding to the geometry of 
a Hall bar will be discussed in section 3. Here the dynamical localization induced by the 
parallel magnetic field will be discussed, as well as the symmetries between the damping and 
the magnetic field. In section 4 the tunneling in the presence of a tilted magnetic field will 
be discussed, with the emphasis on the behaviors near the dynamical localization induced 
by a parallel magnetic field. We briefly summarize in section 5. 



In this section we use the path integral method to reduce the action in the infinite dimension 
to an effective one dimensional one. 

The tunneling is described by the Euclidean action 



A = (0,B±x-B\\z,0) , 



(6) 



2 Effective Action 




dr -mR 2 + i-(B±x — B\\z)y + Viz) + -muj 2 x x 2 + -mLoiy 2 



c 





4 



where (3 = 1/ksT is the inverse temperature. The tunneling rate is equal to exp{— S c /h}, 
where the semiclassical action S c is determined by the bounce solution of the equation 5S = 0, 
in which the periodic boundary condition, (r(h/3), {qj(h/3)}) = (r(0), {^(0)}), is required. 

This is so-called the method of the imaginary part of the free energy F = —ksTlnZ 
calculated from the partition function 



z = J dR'4 q ;.} d(k\ { q ;.} ; r, { q ;» 



with the density matrix 



£>(R'(0), { q ;.(0)};R"(^),{ q ; / (^)}) = /PRP{%} e X p{-~£[R,{ q ,}] 



(9) 



This method is identical to the WKB method at zero temperature, and allows us to have a 
unified treatment of the escape rate for finite temperatures. [§] 

We are interested in the particle tunneling out of the metastable state z = in the 
z-direction. After the tunneling other degrees of freedoms, x and y, as well as {qij}, can 
take any allowed value. Therefore the summation over final states, integrations over the x 
and y, as well as {qij} coordinates, will be taken in the calculation of the partition function. 
Those are Gaussian integrals. To integrate over the x, y and {qij} coordinates, we perform 
a Fourier transformation on the time interval [0, K/3]: 



oo 

(x(T),y(r),z(T)) = — (x n ,y n ,z n )e iUnT 

-oo 



and 



oo 

n l J n=-oo 

Here v n = 27in/h/3 . The action 5* of Eq.(7) can be then rewritten as 



(10) 



(11) 



S[{R n },{qi, j ,n}] = V(z) n + — E 

n P n=-cx 



^mv 2 n z n z„ n + -{mvl + mu 2 x )x n x_ n 



5 



1 e 

-{mvl + mu 2 y )y n y_ n + -( J B_ L x„ - B\\z n )v_ n y_ n 



1 00 

+ ^ ^ ^ 

„=-oo ;j 



1 



2-Ki.n 5iJ,-n - 



with 



/■ft/3 

V(z) n = hp / dry(z(r)) . 
jo 



(12) 



(13) 



Integration over {qij, n }, we have 

1 00 rl 1 

S eff [{R n }] = V(z) n + — o( mZ/ n + iz,n)ZnZ-n + ^{ mV l + muJ l + £x,n)x n X- r . 

">P n=-oo LZ Z 

1 00 rl e 

+ S 9( mZ/ n + m ^+^,«)l/nl/-n + -( 5 -L ;r n- 5 lFn)^-n^ 

n=-oo L2 C 



with 



1 r 00 

&,n = - / duJi(u) 
7T JO 



(14) 



n 1 

Here I — x,y, z. Similarly, the integration over {x n , y n } can be done. The resulting effective 
action is 



rhfi 

efM r )} = / dT 

Jo 



1 1 r h/3 r h/3 

-mz 2 + V (z) + - / dr dr'[k(r - r') + g(r - r')] [z(r) - z(r')] 2 , 
' J 2 Jo Jo 

(15) 

with the normal damping kernel k as 

U \ 1 fj » , coBh[ M (f - |r|)] 

and the damping kernel due to the magnetic field mixed up x — y direction motion contri- 
butions, which we shall call the anomalous damping kernel, as 



1 00 1 



hp ^fr'oo 2 mvl + mw 2 + ^ n 



1 ~ 1 
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[mz/2 + raw 2 + 




i/2 + (ttm/2 + mul + ^,„)(mz/2 + 771^ + £ yn ) 



(17) 
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In obtaining the normal damping kernel k we have used £ n -j^i^ VnT = 

and dropped a term containing the periodic delta function J2^=-oo e%Vnl '■ I n ^ ne large r limit, 

from Eqs.(2) and (16) the normal damping kernel takes the form 



which demonstrates the one-to-one correspondence between the low frequency part of the 
spectral function and the long time behavior of the damping kernel. 

Now we have obtained an effective one-dimensional problem. The physics due to the 
influence of the damping and the magnetic field is buried deeply in the normal and anomalous 
damping kernels. In the following we study two special geometric arrangements to reveal 
the implications. 



In this section we study the tunneling of the charged particle confining to two dimensions, 
that is, there is no perpendicular magnetic field perpendicular to the x — y plane, and the 
^-direction motion of the particle is irrelavent. 

The effective action is given by Eq.(15), with the normal damping kernel of Eq.(16). In 
the absence of the perpendicular magnetic field, B± = 0, the second term in Eq.(17) is zero. 
The anomalous damping kernel g is 




■s+1 ' 



(18) 



3 Edge State Tunneling 




This kernel can be split into two terms, 




h(3 



1 




1 muj 2 y + 

2 mvl + muol + 
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The first term is a delta function, which gives no contribution to the semiclassical action 
under the periodic boundary condition. We only keep the second term: 

air) = -1 ("-EX f 1 (19) 

Thus the anomalous damping kernel takes a particular simple form. This allows us for a 
more detailed calculation of the semiclassical action. 

3.1 Parallel Magnetic Field Only 

This is the case of the tunneling between two parallel equal potential states in the presence 
of a magnetic field, a situation identical to the edge state tunneling in the Hall bar. The 
magnetic field, B\\, is parallel to the x — y plane but perpendicular to the plane of the 
Hall bar. Similar cases have been studied in vortex tunneling in superconductors and and 
superfluids. [7| 

From Eq.(19) in the limit of u y = 0, only the n = mode has a finite contribution. Then 
the effective action is 

S.„ = f d r [- m e + Viz)] + ± (f B|l ) 2 jjjj jf * jf *W - ^)f . (20) 

We note that in the language for the spectral function of the dissipative environment this 
corresponds to the subohmic bath case with s = 0(c.f. Eqs.(2) and (18)), and Eq.(20) is 
explicitly gauge invariant under the change z — > z + constant, because of the periodic 
boundary condition of y imposed in the tunneling calculation. 

For a small enough magnetic field B\\, there is a tunneling solution for the semiclas- 
sical equation SS e ff[q] = 0. Particularly, for a very small B\\ and low temperatures the 
semiclassical action may be evaluated perturbatively: 



S r = I dr 

-oo 



\mil(T) + V(z c (t)) 



+ -i- ( q - B[ 

2m \c 



oo I / roo 



(21) 



where z c {r) is the bounce solution at zero temperature without the magnetic field. Eq.(21) 
shows a remarkable B?> dependence and linear temperature dependence. This feature has 
been verified experimentally. || 

For large enough magnetic fields, the tunneling rate vanishes at zero temperature. One 
way to understand this result is that the magnetic field effectively renormalizes the original 
potential such that the state near z = becomes stable. This can be noted from Eq.(20) 
where as T — > the cross term in the last term vanishes linearly in T, and results in an 
additional term in the potential. The effective action has the following form in this case: 



Thus we have obtained that the s = dissipative environment is marginal for localization 
in tunneling decay, compared to the s — 1 case for the tunneling splitting^]. This is the 
dynamical localization caused by the parallel magnetic field |§. It is worthwhile to point out 
that according to Eq.(23) although a large magnetic inhibits tunneling, a large particle mass 
instead favors the tunneling, contrast to the usual tunneling process where a large particle 
mass disfavors the tunneling. 

An alternative understanding of this dynamical localization is to imagine that the particle 
has a finite cyclotron orbit inside the tunneling barrier. Increasing the magnetic corresponds 
to decreasing the orbit size. Eventually the orbit is smaller than the barrier region, and 
the particle is trapped inside the barrier. Tunneling out of the metastable state becomes 
impossible. However, a finite temperature will destroy this localization. We will return to 
this point in section 4. 




(22) 



A straightforward calculation leads to the criterion for the localization as 




(23) 
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3.2 Finite Dissipation 

Now we set the impurity potential ujy = in Eq.(19). Both the normal damping kernel and 
the anomalous damping kernel are finite. Since the quantum tunneling corresponds to the 
limit of large imaginary time, %f3 — > oo, we look for the large time limit behavior of the 
anomalous damping kernel g. In the this limit, we may replace the summation l/h/3J2 n 
by the integration l/2ir J dv. Then for the environment < s < 2, we find the anomalous 
damping kernel g in the large r limit as 

9 ( T ) = (2 - % ) (£*,)' i-^, (24) 
with a a numerical constant of order of unity. Here 

_ 2 Tjy f 00 A 



^r iz ^_ (25) 

7r m Jo z z + 1 



and Jyiuj) = i] y u Sy has been used. The effective dissipative environment corresponding to 
the anomalous damping kernel is s e ff = 2 — s(c.f. Eqs.(2) and (18)), which leaves the ohmic 
damping unchanged, transforms the subohmic damping into superohmic damping and vice 
versa. 

For the case s > 2, using Eq.(3) for the spectral function J, we find that the effective 
dissipative environment corresponding to the anomalous damping kernel is s e ff — 0, which 
smoothly connects the result for < s < 2. 

An important example is the ohmic damping case. Carrying out a detailed but straight- 
forward calculation, we find the effect action at zero temperature as 

dr -mf + V^y) + — Veff / dr / dr'- -[y(r) - y{r)f , (26) 

-oo 12 J 271 J-oo J~oo \T — T \ 

with the effective damping strength r] e ff as 

*" = <• + (^) 2 4' (27) 
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The effect of a dissipative environment can now be summarized as follows. It has been 
demonstrated in Ref. j| that the subohmic dissipation has strong effects on tunneling, while 
the superohmic dissipation has weak effects. From the above analysis we have that if the 
normal damping kernel k(r) of Eq.(16) is subohmic, the anomalous damping kernel g(r) of 
Eq.(19) is superohmic, and vice versa. Then according to Ref.[[|] the effect of the magnetic 
field, represented by the anomalous damping kernel, is weak/strong on particle tunneling for 
the normal subohmic/superohmic damping. In particular, for normal superohmic damping 
with s > 2, the particle tunnels as if there were no effect of dissipation. For the normal 
ohmic damping with s = 1, we need to compare the relative strength of the magnetic field 
and the dissipation according to Eq.(27). Classically, the Lorentz force tends to keep a vortex 
moving along an equal potential contour, but the friction instead along the gradient of the 
potential. In general, we can conclude that the dissipation tends to suppress the effect of 
the Lorentz force on the tunneling. 

Eqs. (24,25,27) suggest a dual symmetry between the damping and the magnetic field. 
Mathematically, this property arises from treating the y-direction motion as an environment 
for the ^-direction motion, where y and z can be regarded as a pair of canonical variables in 
the presence of the magnetic field, as suggested by the guiding center dynamics, which will 
become more clear in the next subsection. 

3.3 Impurity Potential 

The normal damping kernel vanishes in this situation. The anomalous damping kernel can 
be expressed by hyperbolic functions. We find that the effective action in this case is 



S e ff= I dr 



o 



-mi 2 + V{z) 



+ —( IbA / dr / dr'x 
Am \c J Jo Jo 



^COShK(f -\t-T>\)] _ 

2 sinh[^] 1 
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It is the superohmic case with s e ff = oo, because the effective spectral function has the form 
5(uj — u x ), which has no low frequency mode. Therefore the tunneling rate is nonzero for 
any magnitude of the magnetic field. This result shows that an impurity has a very strong 
influence on the tunneling in the presence of the Lorentz force, because the introducing of the 
impurity potential bends the straight line trajectory of a particle and makes the transition 
to other trajectories possible. We note that by letting u x = we recover Eq.(20), where the 
tunneling rate vanishes for a sufficiently magnetic field. 

We can evaluate the semiclassical action perturbatively, if the magnetic field is weak, 
or, the impurity potential is strong as done in Eq.(21). In the strong magnetic field limit 
semiclassical action may be evaluated by a variational method similar to Ref. Q in the case 
of ohmic damping. However, we have a more powerful method to perform the calculation, 
because a strong magnetic field freezes the kinetic energy of the particle. In this case x and y 
coordinates now form a pair of canonically conjugate variables, the guiding center dynamics, 
and the Euclidean action is|| 



Following the calculation outlined in Ref. 0, we find the semiclassical action in the form 



with Zt the turning point determined by the equation V(y) = V(0). The result shows that as 
mujy increases, the semiclassical action decreases. Therefore we conclude that the impurity 
potential in the ^-direction helps the tunneling in ^-direction. 

3.4 Finite Impurity Potential and Dissipation 

In general we need to go back to Eqs. (15-17) to study the tunneling under the influence of 
an impurity and dissipation. However, based on the insight gained by the above analysis we 




(29) 




(30) 
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can draw a general conclusion: Since both the dissipation and impurity tend to suppress the 
effect of the magnetic field, their total effect will do the same. Indeed the anomalous damping 
kernel is always superohmic. In particular, in the presence of a strong impurity potential and 
ohmic damping, the superohmic-like anomalous damping kernel g may be ignored compared 
to the ohmic normal damping kernel k. The effect of g is to renormalize bare parameters 
such as the mass and the potential. Then from Eq.(15) we may have the effective action as 

S eff = j o dr [-mz 2 + V{z)\ dr jf dr'j-—^[z(r) - z(r')] 2 , (31) 

which looks as if there were no effect of the magnetic field. This result is in accordance 
with the calculation of the change of the tunneling behavior in the presence of random 
impurities. 



4 Near the Dynamical Localization 

In the previous section we have found semiclassically that for a strong enough parallel mag- 
netic field the tunneling rate is zero at zero temperature. This is a dynamical localization 
induced by the magnetic field. Questions remain as regarding to the nature of the transition, 
and the behaviors near this localization. In this section we calculate the tunneling rate near 
the dynamical localization in the absence of the dissipative environment, but allowing a finite 
perpendicular magnetic field and a finite temperature. We will further assume the impurity 
potential is zero to make the analytical analysis transparent. We note that the impurity 
potential behaves as a finite perpendicular magnetic field. 

Without the impurity potential and the dissipation, the effective action is 

S ef f[z(r)}= dr -mz 2 + V(z) +-/ dr dr' 1 
Jo 12 12 Jo Jo 



2 \ mc 



x 



^ coshMf | r -^|)] _ 

2 sinh[^] K ' 

13 



with oo±_ = B±e/mc . The form of this action is identical to the one given by Eq.(28). There- 
fore, the introducing of the perpendicular magnetic field immediately renders the tunneling 
rate finite. To see the nature of the dynamical localization induce by the magnetic field, we 
let u)±_ = 0, and goes back to the renormalized potential at zero temperature( c.f. Eq.(29)): 

V re (z) = V(z) + ^-<^B^z\r). (33) 

The semiclassical action is always finite as long as the metastable exists. This implies that 
the dynamical localization is a 'first order' phase transition as varying the magnetic field: 
below the critical field B c the semiclassical action is finite, and above B c it is infinite( or 
a jump to another value if the renormalized potential has the stable point not at z — 0). 
Hence the first set of result is 



S c 



S c0 T = 0, B ± = 0, B\\ < B c 

(34) 

oo T = 0, B ± = 0, B\\ > B c 



Here S c o is the semiclassical action right at the transition. 

For a finite but small temperature and zero perpendicular magnetic field, the effective 
action is 

r1 1 m ( eB»\ 2 1 r h P r h P 

If B\\ < B c , the temperature can be treated as a perturbation. The calculation is identical 
to the one leading to Eq.(21). If B\\ > B c , the temperature dependence is singular. The 
semiclassical action may be estimated in the following way. In this case the renormalized 
potential has two local minima, that is, V'. e (z) = has two solutions (We only consider such 
a situation here. For a more complicated potential V, more than two minima is possible.), 
one at z = and one at zi determined by the magnetic field. The second minimum is 
metastable. The tunneling process can now be viewed as two steps. First, by the thermal 

14 



activation the particle raises its energy to that of the second minimum. Then it tunnels into 
the second minimum. Therefore the semiclassical action has the form of S c0 + (V re (z 2 ) — 
V re (0))/kBT. Here we note (V re (z2) — V re (0))\B n =B c = 0. Near the dynamical transition 
point, the semiclassical action can be summarized as 



(35) 



S c o B± = 0, B\\ < B c 

h{Vre{ f- T Vrem + S c0 B ± = 0, B\\ > B c 
There is an Ahrenius-type temperature dependence, with a small effective activation energy 
V re (z 2 )-V re (0)- 

In the case of finite perpendicular magnetic field, B± ^ 0, the general expression should 
follow from Eq.(32). Near the dynamical localization point, again simple expressions for the 
semiclassical action can be obtained. If the temperature is high such as ksT > hu±/2, the 
effect of B± can be ignored. In this case the semiclassical action takes the same form as for 
B± = 0, and is 



S c = 



Sco 2 A; B T > hu ± , B\\ < B c 

h(VM-v re (o)) + ScQ 2k B T>hw ± ,B ll >B c 

In the opposite limit, k B T < hu±/2, the effective action may take the following form: 



(36) 



rH/3 

[z(r)\ = I • 
Jo 



dr 



-mz 2 + V{z) 



m ( eB\\ 
4 \ mc 



j^dr dT'^e~^-^[z(r)-z(r')f 



If B\\ < B c , the semiclassical action can be evaluated perturbatively. For B\\ > B c , we make 
following estimation. Since fi(3 > 2/cj_l, u± as a cutoff may be introduced. The effective 
action becomes: 



'e// 



Jo 



dr 



-mz 2 + V(z) 



f d \ 2 2 2 



+ 4 



mc 



dr dr'^[z(r) - z{r')f . (37) 



This suggests that u± plays the same role as the temperature, and the semiclassical action 
will have the form S c0 + (V re (z 2 ) — V re (0))/tUv±. We summarize the results in this low 
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temperature limit, 

f S c0 2k B T < hw x ,B\\ < B c 

S c ={ . (38) 

[ 2(Ke(22 J~ Ke(0)) + S c0 2k B T < huj ± , B\\ > B c 

In the small uj± limit S c goes to oo when B» > B c . 

We must point out that the calculations in this section are of mean field type. No 
fluctuations around the mean field solutions are considered, which may be significant in 
some situations. 



5 Summary 

We have performed a systematic study of the influence of the dissipation and the magnetic 
field on tunneling in the presence of an impurity potential. Along the tunneling coordinate 
the effective action is one dimensional, but with normal and anomalous damping kernels 
which are no-local in time. Two cases, the Hall bar geometry and the behaviors near the 
dynamical localization, have been investigated in detail in order to reveal the physics implied 
in the effective action. A dual symmetry between the damping and the magnetic field has 
been found, corresponding the mapping between different heat baths. The very pronounced 
result is the dynamical localization caused by the parallel magnetic field: For a strong enough 
parallel magnetic field the tunneling can be completely suppressed. For the weak parallel 
magnetic field magnetic the tunneling rate is finite, where the semiclassical action quadrati- 
cally depends on the strength of the magnetic field and linearly on the temperature. Explicit 
expressions in various limits have been obtained, which can be directly tested experimentally. 
For example, for the system of electrons on a helium surface^, the critical magnetic field for 
the dynamical localization is about 1 Tesla, and for the heteroj unction systems it is about 
50 Tesla, depending on the material parameters. 
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